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On the nature of the circumstellar medium of the 
remarkable type la/IIn supernova, SN 2002 ic 
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ABSTRACT 

We present results from the first high resolution, high S/N, spectrum of SN 2002 ic. 
The resolved Ha line has a P Cygni-type profile, clearly demonstrating the presence 
of a dense, slow-moving (~100kms~^) outflow. We have additionally found a huge 
near-IR excess, hitherto unseen in type la SNe. We argue that this is due to an 
IR light-echo arising from the pre-existing dusty circumstellar medium. We deduce 
a CSM mass probably exceeding 0.3Mq produced by a mass loss rate greater than 
several times 10~^ MQyr~^. For the progenitor, we favour a single degenerate system 
where the companion is a post-AGB star. As a by-product of our optical data, we are 
able to provide a firm identification of the host galaxy of SN 2002 ic. 

Key words: circumstellar matter - supernovae: general - supernovae: individual: SN 
2002 ic stars: winds, outflows - dust 



1 INTRODUCTION 

Tlie lack of a convincing detection of hydrogen in the 
spectrum of any type la supernova has been difficult to 
reconcile in the single-degenerate scenario where the non- 
degenerate companion, in most candidate progenitor chan- 
nels, is hydrogen- rich. H[owever, Hamuy et al. (2003) recently 
announced the discovery of }ia emission associated with the 
type la supernova, SN 2002 ic. Over a time-span of about 
+7 to +48 d from maximum light (tsmax = JD 2452601 = 
2002 November 22; Hamuy et al. 2003, Fig. 4), the optical 
spectra of SN 2002 ic exhibited similar but weaker features 
to those of 'normal' type la SNe. However, strong Ha emis- 
sion was also apparent: the Hq feature consisted of a narrow 
component (unresolved at 300kms~^) atop a broad compo- 
nent (FWHM ~1800kms~^). While the narrow component 
could have beeen due to an underlying HH region (but see 
below), Hamuy et al. (2003) argued that the broad compo- 
nent arose from ejecta/circumstellar medium (CSM) inter- 
action, and that this interaction also provided the continuum 
source required to dilute the spectral features of SN 2002 ic. 
By day -1-48, they found that the spectrum could be equally 
well-matched by either a suitably 'diluted' coeval spectrum 
of the type la SN 1990 N, or by an unmodified, roughly^ 
coeval spectrum of the type Tin SN 1997 cy. Type Hn su- 
pernovae (SNe Hn) are so called because of the presence at 
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^ Note that the explosion date of SN 1997 cy is uncertain Ger- 
many et al. (2000). 



early times of narrow lines in the spectra originating in a rel- 
atively undisturbed circumstellar medium (CSM) (Schlegel 
1990). Their progenitors must therefore have undergone one 
or more mass-loss phases before explosion. 

In order to investigate the origin of the hydrogen emis- 
sion and hence the nature of SN 2002 ic and its circumstel- 
lar environment, we have acquired high resolution optical 
spectroscopy at -1-256 days, and HK-hand IR photometry 
at -1-278 and 380 days. The first results of this study are 
presented here. 



2 OBSERVATIONS 

2.1 Optical Spectroscopy 

We obtained optical spectra of SN 2002 ic and its pur- 
ported host galaxy on 2003 August 05 (— +256 d) with the 
ESQ Very Large Telescope (VLT) Unit 2 (Kueyen) and 
Ultra- Violet Echelle Spectrograph (UVES). We used a 3" 
(PA=90°) slit which yielded a resolution of ~9kms~^. The 
seeing was ~ 0'.'9. The exposure times for SN 2002 ic and the 
galaxy were 2200 s and 1100 s respectively. The data were re- 
duced in the Figaro 4 environment. Wavelength calibration 
was by means of a ThAr arc taken at the end of the expo- 
sure of each of the targets. Flux calibration was carried out 
with respect to the spectrophotometric standard Feige 110. 

A portion of the UVES spectrum obtained on day 256 
is shown in Fig. 2.1. The spectrum is dominated by a strong 
Ha feature. In addition, a weak broad feature around 9000 A 
is present which is the blend of O I 8446A and the Ca II IR 
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Figure 1. (a): A portion of tiie pipeline-reduced red spectrum of SN 2002 ic at +256 d in tlie restframe of tlie supernova. Tlie gap at 
8600 A is due to a gap between the mosaic of the two red chips. Above Fig. la is shown a portion of the pre-extraction 2-D spectrum. 
Ha emission from galaxy A is marked by 'A' (see also Fig. 3 and Sec. 3); the supernova continuum is marked by 'SN'. (b): Expanded 
view of the unbinned Ha profile, (c): Evolution of the Ha profile of SN 2002 ic. The numbers on the left-hand side give the epoch (days) 
while the numbers on the right (in brackets) give the factors by which the profiles have been shifted vertically for clarity. The top 4 
spectra are from Hamuy et al. (2003); the 217 d spectrum is from Deng et al. (2004); the 256 d spectra (thick red lines) are derived from 
our UVES spectrum, binned to match the resolutions of the Hamuy et al. (2003) and Deng et al. (2004) spectra repectively. 



triplet (Wang et al. 2004). In Fig. 2. lb we show the Ha profile 
in more detail: it comprises a narrow, but resolved, P Cygni- 
like profile atop a broad emission feature. There may also be 
a very broad feature present, but owing to the blue limit of 
the spectrum, we are unable to give a complete description of 
this feature. The Deng et al. (2004) -1-217 d spectrum extends 
further to the blue, and they attribute the very broad fea- 
ture to [O I] 6300,6364 A, (FWHM ~ 26000 kms"^). We do 
not give further consideration to this component. However, 
we note that as a consequence of the presence of the very 
broad feature, some authors refer to the "broad emission 
feature" mentioned above and shown in Fig. 2.1 as the "in- 
termediate component" . We shall continue to refer to this as 
the "broad" feature. It is about 5800 km s~^ across the base 
(FWHM~1550kms^^). In order to extract more detailed in- 
formation about the narrow feature, we generated a model 
P Cygni profile using a homologously expanding CSM above 
a photosphere, with a rest frame wavelength equal to that 
of the narrow peak and an exponentially-declining density 
profile. The P Cygni profile parameters were adjusted by 
eye to match the absorption component, yielding a velocity 
of 100 km s^^ at the photosphere, and an e-folding velocity 
of 30kms~\ The maximum detectable extent of the blue 
wing of the absorption is ~250kms~^. The P Cygni model 
also demonstrated that the narrow component includes ad- 
ditional emission not taken into account by the absorption. 
Both the narrow emission component and the P Cygni pro- 
file suggest a CSM velocity of 80-100 km s~^. In addition, 
there is a small but significant shift between the narrow 
emission component and the 1500 km s~^ component in the 
sense that the narrow component peak is 111 ±7 km s~^ fur- 
ther to the red. The Ha profile parameters are summarised 
in Table 1. 

In Fig. 2.1c we show the evolution of the Ha profile 
from +7 to +256 d. Our +256 day UVES spectrum provides 
a good match to the +217 d profile, showing that the narrow 



Table 1. Parameters for the Ha profile of SN 2002ic at +256 d. 



Narrow Broad 
Peak wavelength (A) 6998.10+0.05 6995.5+0.15 
Width (FWHM A) 1.83+0.04 36+0.5 
Width (FWHM km/s) 78.4+1.7 1500+20 
Intensity (10~^^ ergcm'^s'^) 4.9+0.1 50.7+0.5 



absorption feature which we detected would not have been 
apparent in the Deng et al. (2004) spectrum. In contrast, 
however, our +256 d UVES spectrum binned to the Hamuy 
et al. (2003) resolution still shows evidence of the narrow 
absorption feature, whereas no such feature is evident up 
to +48 d. Indeed, there is little evidence of an asymmetry in 
the profile at early times. This suggests that the photosphere 
is located in a thin, opaque, cool dense shell (CDS) formed 
by the reverse shock propagating into the extended stellar 
envelope (e.g. Chevalier & Fransson 1994). 

2.2 Near-infrared photometry 

We obtained a /C-band image of SN 2002 ic at UKIRT on 
2003 August 27 (= +278 d) using UIST; subsequent H and 
if-band images were obtained on 2003 December 7 (+380 d) 
using UFTI. The data were reduced using standard proce- 
dures implemented in the pipeline software. The 278 d K- 
band image is shown in Fig. 3. The supernova is clearly vis- 
ible 5" west of galaxy A. 

Magnitudes were measured using aperture photometry 
within the Starlink package GAIA. The sky background was 
measured using a concentric annular aperture. An aperture 
radius of 0'.'72 was selected, equivalent to 6 and 8 pixels for 
UIST and UFTI respectively. The sky annulus was chosen to 
have inner and outer radii, respectively, xl.5 and x2.5 that 
of the aperture. Magnitudes were determined by comparison 
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Figure 2. Decomposition of the Ha profile of SN2002ic at 
+256 d. The thick line shows the combination of the three Gaus- 
sians and the P Cygni model (dashed-dotted line). 

Table 2. Infrared magnitudes of SN 2002ic 



Date day H K H - K 
20030827 +278 ■ ■ ■ +18.01+0.03 

20031207 +380 +18.98+0.02 +17.76+0.03 +1.22+0.04 
Notes. The errors shown are statistical only. 

with the UKIRT standard star FS 105. We assess the overall 
uncertainty in the photometry to be about ±5%. The HK 
magnitudes are listed in Table 2. 



3 REDSHIFT OF SN 2002IC AND HOST 
GALAXY IDENTITY 

We determined the redshift of SN 2002 ic from the narrow 
emission component of the P Cygni profile and found it to 
be 2 = 0.0663. This is consistent with the measurement 
of Hamuy et al. (2003). Our redshift implies a distance of 
~280 Mpc {Ho = 70kms"^ Mpc"^). According to Hamuy 
et al. (2003), the redshift of the galaxy ~ 5" E of SN 2002 ic 
(marked A in Fig. 3) is 0.22 thereby ruling out association 
with the supernova. The only other nearby galaxy is the one 
marked B in Fig. 3, lying ~ 10" S of the supernova; Hamuy 
et al. (2003) suggested an association but do not report a 
redshift for this galaxy. However, our UVES spectrum of 
galaxy B indicates z = 0.0784, i.e. it is unlikely to be the host 
galaxy. Fortuitously, during our UVES observation of the 
supernova, the extreme eastern end of the slit intercepted 
the nuclear region of galaxy A. We noticed an emission fea- 
ture at the corresponding spatial position (along the slit) in 
the spectra, and in the same order as the Ho? feature of SN 
2002 ic, shifted only slightly in wavelength. Assuming the 
feature was also due to Ha emission (but from Galaxy A), 
we derive z — 0.0667. This indicates that Galaxy A must be 
the host of SN 2002 ic. Hamuy (priv. comm.) confirms that, 
owing to a target acquisition error, the redshift given for the 
host galaxy in Hamuy et al. (2003) is incorrect. 
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Figure 3. X-band image of SN 2002ic obtained with UIST at 
UKIRT on 2003 August 27 (+278 d). The scales give the off- 
sets relative to SN 2002 ic: aj2000 = 0l''30'"02^55; i5j2ooo = 
+21°53'06".9. We argue in Sect. 3 that the galaxy marked 'A' is 
the host galaxy. 

4 DISCUSSION 

4.1 Constraints from the Ha profile 

The narrow component of the Ha feature suggests an ori- 
gin in a wind fiowing at Vm ~ lOOkms"^. The presence and 
velocity of the P Cygni-like absorption immediately rules 
out an origin in a line-of-sight HII region i.e. the narrow 
emission/absorption feature is intrinsic to the supernova or 
its immediate environment. The similarity of the early-time, 
low-resolution spectra of SN 2002 ic to that of the type Iln 
SN 1997 cy has been noted by several authors. We find that 
the similarity also holds at high resolution; the late-time 
Ha profile of SN 2002 ic is compared to those of type Tin su- 
pernovae observed at comparable epochs in Fig. 4. Narrow 
emission/absorption profiles superimposed on broad emis- 
sion features have been observed in the type Tin events SN 
1997 ab at 425 d and SN 1997 eg at +202 d (Salamanca et 
al. 1998, 2002). For SN 1997 ab the narrow absorption blue 
wing limit yields a velocity of 90kms~^, superimposed on 
an emission feature of ~1800kms~^ FWHM, with wings 
extending to ~4000kms~^. For SN 1997 eg, the correspond- 
ing figures are 160kms~^, 3800kms~^and ~11000kms~^. 
In both cases the narrow feature is displaced redward of 
the peak of the broad component by ~600kms~^. Sala- 
manca et al. (1998) attribute the apparent relative shifts to 
self-absorption in the intermediate component which pref- 
erentially attenuates the red wing. As indicated above, the 
SN 2002 ic narrow component also exhibits a 'redshift' rela- 
tive to the broad peak, but at lllkms"^ the shift is much 
smaller, suggesting that self-absorption is less important. 
The luminosity of the narrow and intermediate Ha fea- 
tures of SN2002ic are, respectively, 5 x lO^^ergs"^ and 

5 X 10*° ergs-^ 

The 1500 km s~^ component is probably produced by 
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Figure 4. Compajrison of the Ha profile of SN 2002 ic with high- 
resolution spectra of two type Iln supernovae at roughly compa- 
rable epochs. The data for SN 1997 ab and SN 1997cg arc taken 
from Salamanca et al. (1998, 2002) and have been corrected for 
their respective redshifts. 



the supernova ejecta/wind interaction, as suggested by 
Hamuy et al. (2003). We note that the width of the broad 
component declined from 1800 km to 1500 km s~^ be- 
tween days +47 and -1-256. This gives additional weight to 
the ejecta/wind interaction scenario. 

The high latc-timc luminosity of SN 2002ic allows us to 
rule out radioactive decay of ^®Ni as the dominant energy 
source. The UBVRI light curve Deng et al. (2004) gives a 
luminosity of 5 x 10*^ crgs~^ at 250 d, whereas the total 
radioactive luminosity of O.TM© ''^Ni at this epoch is only 
1 X 10*2 ergs"^ Moreover, by +250, only about ~10% of 
the decay gamma-rays would be deposited in the ejecta of 
the presumed type la SN ejecta (e.g. Axelrod 1980). There- 
fore, the dominant source of energy for the broad component 
luminosity, L^°^'^, must be due to the ojccta/CSM interac- 
tion. In this scenario, L^°'^'^ is proportional to the kinetic 
energy dissipation rate across the shock front. We can use 
this luminosity to estimate the mass-loss rate, M, (Sala- 
manca et al. 1998): 



L Broad 1 
Ha — -l^H. 



M 



(1) 



where is an efficiency factor which peaks at ~0.1. Vs 
is the shock velocity, and Vw is the velocity of the un- 
shocked wind, assumed to be freely-expanding. Prom the 
broad Hq line, Va= 2900 kms~^, while the narrow feature 
gives tJ™ = 100 kms^^. Substituting into the above equation, 
and using the broad component luminosity, wo obtain a 
mass-loss rate of M ~ 1.2 x 10"2/(eH„/0.1) Moyr"^ 

Salamanca et al. (1998) also show that by using the 
luminosities of both the broad and narrow omission com- 
ponents, the CSM density and mass can be estimated. Us- 
ing the component luminosities given above, we find that 
the number density at the CSM inner limit on day 256 is 
Ui = 1.35 X lO'^cm"^. Assuming that the CSM was cre- 
ated by a steady wind (density oc r~^) and that the in- 
ner limit of the CSM on day 256 corresponds to the radius 
reached by the 2900 km s~^ shock, we find a total CSM mass 
of O.O37iZ2/-RiM0, where Ri and R2 are, respectively, the in- 



ner and outer limits of the CSM, and R2 » Ri- Ri can be 
identified with the shock radius. 



4.2 Constraints from NIR photometry 

We now consider the IR emission. On day +380, H — K = 
1.22±0.04. Such a colour corresponds to a 1430±40K black- 
body. We therefore propose that the late-time IR emission 
is due to thermal emission from hot dust associated with 
SN 2002 ic or its progenitor. However, it is possible that 
the IR flux contains a component due to hot (T~10,000K) 
residual photospheric emission such as might be produced 
by the sliock/CSM interaction. We therefore measured the 
continuum level in the vicinity of 0.94/im on day +256 and 
extrapolated to day +380 assuming an exponential decline 
timescale of 170 d (Deng et al. 2004). We then extrapolated 
to the H and K bands assuming a Rayleigh- Jeans law. Prom 
this we conclude that ~30% and ^^^8.5% of the H and K 
band fluxes respectively was due to contamination by the 
hot photosphere. After subtracting the photospheric compo- 
nent and correcting for a Galactic extinction of Ay = 0.198 
(NED), we find that the net IR fiux on day 380 can be 
reproduced by a T=1220K blackbody and a luminosity of 
The corresponding figure for +278 d (as- 



1 X lO'*^ ergs"^ 



suming the same temperature) is 0.7 X 10 ergs" . We note 
that if the IR emission were due to dust condensation in 
the ejecta, for the corresponding 1220 K blackbody to at- 
tain these luminosities it would need to have expanded at 
8000-9000 km s"'^ since the supernova exploded. 

We now consider the location and origin of the hot dust. 
We first note that as with the UBVRI emission, contempo- 
rary radioactive decay can have made only a minor contri- 
bution to the late-time IR luminosity of SN 2002 ic. We also 
note that dust is not expected to condense in type la explo- 
sions. This is consistent with the fact that to produce the 
observed IR luminosity, the 1220 K blackbody surface would 
have to be located as far out as the ~8500kms~^ region of 
the ejecta. Dust condensation in such circumstances seems 
unlikely. We conclude that the IR luminosity arises from a 
pre-existing dusty CSM. Heating of this dust can be via (a) 
local heating by the ongoing ejecta/CSM shock interaction, 
or (b) photon emission from the supernova yielding an IR 
echo in the unshocked CSM. Deng et al. (2004) find veloc- 
ities exceeding 10000 km s"'^ in C, O, and Ca, which they 
attribute to the supernova ejecta. Thus, local heating of a 
dusty CSM by the ejecta/CSM shock might seem to be a 
possibility. However, the the initial fiash from the supernova 
would have evaporated any CSM dust to at least ~3500 AU 
for carbon-rich grains (Tovap ~ 1900 K) and 16,500 AU for 
oxygen- rich grains (Tcvap = 1500 K) (Dwek 1983, 1985). 
Yet a 10000 km s^^ shock would have reached only 1600 AU 
by +278 d. We therefore rule out local heating of the dust 
by the ejecta/CSM shock. 

We now tost the possibility that the IR omission arose 
from the heating of CSM dust by photons from the super- 
nova i.e. the IR-echo scenario. We use the bolometric light 
curve of Deng et al. (2004) which can be approximately de- 
scribed as having a peak (t=0) luminosity of Lo = 3 x 10*^ 
ergs^^ and an exponential decline timescale of 170 d. We 
attribute this slow decline to the energy released in the in- 
teraction of the shock with the CSM. The bolometric light 
curve is based on UBVRI photometry, and so does not in- 
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elude possible additional energy from an X-ray precursor. 
However, given the much higher opacity of dust grains to 
UV-optical light, it is likely that the X-ray contribution to 
grain heating will be small. 

Dwek (1983) showed that the IR-ccho light curve com- 
prises an initial plateau phase, followed by a decline. The 
transition from plateau to decline corresponds to the passing 
of the ellipsoid vertex from the dust-free cavity into the re- 
gion of unevaporated dust. Thus, the radius of the dust free 
cavity Ri = cti /2. We noted that the IR flux from SN2002 ic 
barely changed between days 278 and 380. From this we con- 
clude that the vertex was still within the dust-free cavity on 
day 380, implying a cavity radius of _Ri > 32800 AU. Wc use 
Eq. 17 in Dwek (1983) to estimate the optical depth of dust 
required to yield the observed IR flux from SN 2002 ic. The 
parameters adopted were: Ri = 32800 AU, d = 280 Mpc, 
dust IR emissivity proportional to for wavelengths down 
to 0.2 /im, a mean UV- visual absorption efiicicncy of 1, and 
an initial dust temperature at the cavity boundary assumed 
to be about equal to an evaporation temperature of 1500 K. 
The CSM was assumed to have been produced by a steady 
wind so that the density is proportional to r~^. We find that 
the IR flux at both epochs is reproduced with a dust opti- 
cal depth of Td ~ 0.01. For a gas-to-dust mass ratio of 160, 
a grain material density of 3gcm~* and a grain radius of 
0.1 /ttm, this translates (from Dwek 1983) into a total CSM 
mass (including the dust-free cavity) of ~0.3J?2/J?i M©, 
where R2 is the outer limit of the CSM and R2 » Ri- Thus 
the CSM mass exceeds 0.3 Mq. The corresponding mass loss 
rate (again for R2 » Ri), assuming a wind velocity of 
lOOkms-i is 1.9 x 10"'' Moyr"^ (Dwek 1983, Eq. 19). It 
is about x5 the value which Dwek (1983) derived for the 
type IIL SN 1979 C. 

Compared with the derivations from the Ha line, the 
IR analysis produces a lower mass-loss rate (~1%). This 
discrepancy could, in part, be due to an underestimate of 
the shock velocity leading to an overestimate of the mass- 
loss rate derived from the broad Ha line. Nevertheless, both 
analyses indicate a CSM mass probably exceeding 0.3Mo 
produced by a mass loss rate greater than several times 
lO"*M0yr-^ The mass-loss rate inferred from this work 
and that of others is higher than expected from traditional 
mass loss mechanisms. We remark that the high values of 
M are - at least partly - a consequence of simplifying as- 
sumptions e.g. clumped winds would mimic a high mass-loss 
rate. 

The close similarity between SN 2002 ic and type Iln 
SNe, has raised doubts as to whether SN 2002 ic is a bona fide 
la event and whether other type Iln (i.e. core-collapse) SNe, 
only discovered at late epochs, may have been SN 2002 ic- 
like events. We suggest that the type Iln phenomenon is 
predominantly related to the amount of CSM around the 
progenitor system rather than the type of explosion. 

A type 1.5 scenario i.e. the explosion of a single, mas- 
sive AGB star has been invoked as a possible progenitor of 
SN 2002 ic (e.g. Hamuy et al. 2003). We point out that an 
extremely low metallicity would be required to inhibit mass- 
loss so as to allow the degenerate core to grow to Mch e.g., 
a 4M0 star would need to have logZ/ZQ < —3 (Ziljstra 
2003). Also, a low mass-loss rate is at odds with the large 
amount of CSM inferred for this event, although we can- 
not rule out the possibility that the CSM is due to a binary 



companion. Furthermore, SN 2002 ic exhibited an exception- 
ally high maximum V-band luminosity and a much slower 
post-~-f-25days decline rate in BVI than is seen in normal 
SNe la. The opposite effect would be expected for type 1.5 
events i.e. the photometric and spectral evolution should be 
similar to II events at early times and dominated by the de- 
cay of radioactive Ni and Co at late times as for type I and 
IIP events (Iben & Renzini 1983). 

Any progenitor scenario must satisfy all the observa- 
tional constraints viz: the type la-likc behaviour at early- 
times and the type Iln behaviour at late-times (Hamuy et 
al. 2003); broad profiles of Ca and O (Deng et al. 2004); an 
asphcrical CSM (Wang et al. 2004); a slow- moving outfiow 
at ~ 100 km s^'^ and a dusty CSM. It must also explain the 
apparent rarity of SN 2002 ic-like events. Taking these obser- 
vational constraints at face- value, we currently favour a sys- 
tem involving a post- AGB star. There are several known ex- 
amples of post- AGB stars that have high inferred mass-loss 
rates and dusty discs CSM) (e.g. IRAS 08544-4431, Maas 
et al. 2003). These objects have typical outflow velocities of 
the order of 100 km s"'^. Furthermore, the post- AGB phase 
can be relatively short, ~10^ — 10* yrs. (van der Veen et al. 
1994). Further planned observations will no doubt provide 
more clues as to the previous evolution of SN 2002 ic. 
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